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ADDITIVE  SCHWARZ  METHODS  FOR  ELLIPTIC 
FINITE  ELEMENT  PROBLEMS  IN  THREE  DIMENSIONS 

Maksymilian  Dryja'  and  Olof  B.  Widlund^ 

Abstract.  Many  domain  decomposition  algorithms  and  certain  multigrid  methods  can 
be  described  and  analyzed  as  additive  Schwarz  methods.  When  designing  and  analyzing 
domain  decomposition  methods,  we  encounter  special  difficulties  in  the  case  of  three  di- 
mensions and  if  the  coefficients  are  discontinuous  and  vary  over  a  large  range.  In  this 
paper,  we  first  introduce  a  general  framework  for  Schwarz  methods.  Three  classes  of  appli- 
cations are  then  considered:  certain  wire  basket  based  iterative  substructuring  methods, 
Neumann-Neumann  algorithms  with  low  dimensional,  global  subspaces  and  a  modified 
form  of  a  multilevel  algorithm  introduced  by  Bramble,  Pasciak  and  Xu. 
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1.  Introduction.  In  this  paper,  we  discuss  additive  Schwarz  methods  for  solving 
systems  of  hnear  algebraic  equations  which  result  from  finite  element  approximations  of 
second  order,  elliptic  problems  in  three  dimensional,  bounded  regions.  A  general  frame- 
work is  presented  which  is  quite  useful  in  the  design  and  analysis  of  a  variety  of  domain 
decomposition  and  some  multigrid  methods.  Three  methods  are  then  described  and  ana- 
lyzed. They  are  extensions  of  methods  previously  considered  in  the  literature  for  solving 
the  systems  of  algebraic  equations  which  correspond  to  the  interfaces  between  the  substruc- 
tures. We  also  consider  problems  with  discontinuous  coefficients  with  a  great  variation  in 
the  values. 

The  first  method  considered  is  an  iterative  substructuring  algorithm  recently  devel- 
oped by  B.  Smith  [9],  [10].  The  second  is  a  domain  decomposition  method  developed  by 
R.  Glowinski,  P.  Le  Tallec,  Y.-H.  De  Roeck  et  al.,  see  [1],  [3].  Finally,  we  consider  a  variant 
of  a  multigrid-like  method  discovered  by  J.  H.  Bramble,  J.  E.  Pasciak  and  J.  Xu  [2]. 

The  paper  is  organized  as  follows.  In  Section  2  a  finite  element  approximation  of  second 
order,  elliptic  problems  is  considered.  A  system  of  algebraic  equations  corresponding  to 
the  discrete  problem  is  reduced  to  a  system  defined  on  the  interfaces  of  the  substructures 
after  eliminating  the  interior  variables  associated  with  the  interior  nodal  points  of  each 
substructure.  This  reduced  sj^stem  with  the  Schur  complement  matrix  can  be  solved  by  a 
variety  of  iterative  methods. 

In  Section  3,  a  general  abstract  framework  for  additive  Schwarz  methods  is  introduced, 
see  also  [5],  [6].  In  Sections  4,  5  and  6  the  three  algorithms  mentioned  above  are  described 
and  analyzed  inside  this  framework.  In  pai'ticular,  we  discuss  problems  with  discontinuous 
coefficients  and  show  that  the  rate  of  convergence  of  certain  variants  of  the  methods  can 
be  made  independent  of  the  variation  of  the  coefficients. 

The  results  of  Section  4  of  this  paper  have  been  obtained  jointly  with  Barry  Smiths; 
see  further  [4]. 

2.  Differential  and  finite  element  model  problems.  To  simplify  the  presentation, 
we  discuss  only  two  model  problems,  a  standard  Poisson  equation  and  a  special  second 
order  problem  with  discontinuous,  piecewise  constant  coefficients.  We  call  them  model 
problem  I  and  II,  respectively.  The  continuous  model  problem  I  is  of  the  form: 

Find  u  e  H](^)  such  that 

(2.1)  a{u,v)  =  f{v),  v^H'oiQ) 
where 

(2.2)  a{u,v)=  I  Vu  ■  Vvch  ,      f(v)  =   f  fvdx 

JQ  ■  Jn 

For  simplicity,  we  let  Q  be  a  bounded  polyhedral  region  in  three  dimensions.  A  coarse 
triangulation  of  Q.  is  introduced  by  dividing  the  region  into  nonoverlapping  simplices  fi,, 
i  =  1, . . . ,  A'',  which  are  also  called  substructures.  The  substructures  fi,  are  further  divided 
into  elements  Cj.  We  associate  parameters  H  and  /;  with  the  coarse  and  fine  triangulations 
and  assume  that  these  triangulations  are  shape  regular  in  the  sense  common  to  finite 
element  theory. 


Let  V'^(Q)  be  the  finite  element  space  of  continuous,  piecewise  linear  functions  defined 
on  the  fine  triangulation  and  which  vanish  on  dQ,  the  boundary  of  Q.  The  discrete  model 
problem  I  is  of  the  form: 

Find  Uh  G  V'^  such  that 

(2.3)  a{uH,Vh)  ^  f(vk)  ,  vhEV' 
or  alternatively,  find  the  vector  x  of  nodal  values  such  that 

(2.4)  A' a-  =  6 

Here  A'  is  the  stiffness  matrix,  x  is  the  vector  of  nodal  values  and  b  the  load  vector. 

Our  goal  is  to  describe  and  analyze  iterative  methods  for  solving  (2.2)  which  can  be 
regarded  as  additive  Schwarz  methods. 

All  our  results  can  be  extended  to  general  conforming  finite  element  approximations 
of  any  self-adjoint,  second  order,  elliptic  problem.  This  includes  the  case  when  there  is  a 
great  variation  in  the  values  of  the  coefficients.  Here  we  consider  only  the  case  when 

(2.5)  a(u,v]==   f  p(x)Vu-Vvdx 

Jn 

where  p{x)  >  0  is  a  piecewise  constant  function.  Equation  (2.1)  with  the  bilinear  form  (2.5) 
is  called  model  problem  H.  We  assume  that  the  jumps  of  p{x)  occur  only  at  substructure 
boundaries.  Thus,  p{x)  =  pi  =  constant  >  0  on  the  substructure  fi,.  The  methods 
discussed  can  be  generalized  to  the  case  when  p{x)  varies  continuously  in  each  subregion. 

Let  K  be  the  stiffness  matrix  given  by  the  bilinear  form  (2.2)  or  (2.5).  In  the  first  step  of 
many  iterative  substructuring  methods  the  unknowns  in  the  interior  of  the  substructures 
are  eliminated.  This  reduces  the  system  (2.4)  to  a  system  of  linear  algebraic  equations 
associated  with  the  interfaces  only.  We  now  describe  this  procedure. 

Let  A'^''  be  the  stiffness  matrix  of  the  bilinear  an,(uft,u/i)  which  represents  the  contri- 
bution of  the  substructure  Q,  to  the  integral  an(u/,,u/i)  =  a{uh,Vh)-  Let  x  and  y  be  the 
vectors  of  nodal  values  that  correspond  to  the  finite  element  functions  Uh  and  v^,  respec- 
tively. Then  the  stiffness  matrix  A'  of  the  entire  problem  can  be  obtained  by  using  the 
method  of  subassembly, 


(2.6)  x^Ky  =  Yl 


M)' 


Here  x^''  is  the  subvector  of  nodal  parameters  associated  with  Q,,,  the  closure  of  Q;.   We 
represent  AT*''  as 

;-(.)         jA') 


(2.7)  1   ^.{,It 


IB 

(«) 

IB        ^'^BB 


dividing  the  subvector  x*''  into  two,  Xj  and  Xg  ,  corresponding  to  the  variables  which  are 
interior  to  the  substructure  and  those  which  are  shared  with  other  substructure,  i.e.  they 
are  associated  with  the  nodal  points  of  5^,.    Since  the  interior  variables  are  associated 


with  only  one  of  the  substructures,  they  can  be  ehminated  locally  and  in  parallel.  The 
resulting  reduced  matrix  is  a  Schur  complement  and  is  of  the  form 

(2.8)  S^'^^K^\,-K\fK\'ri^\i 

From  this  it  follows  that  the  Schur  conrplement  corresponding  to  the  global  stiffness  matrix 
K  is  given  by  S  where 

(2.9)  xlSys  =  Y.^B^S^'yB^ 

If  the  local  problems  are  solved  exactly  what  remains  is  to  find  a  sufficiently  accurate 
approximation  of  solution  of  the  linear  system 

(2.10)  SXB  =  98 

Note  that  the  elimination  of  the  interior  variables  of  a  substructure  can  be  viewed  in  terms 
of  an  orthogonal  projection,  with  respect  to  the  bilinear  form,  of  the  solution  Uh  of  (2.3) 
onto  the  subspace  i?o(^t)  ^  ^'^''-  ^^  ^^  easy  to  show  that  these  subspaces  are  orthogonal,  in 
the  sense  of  a(u,v)^  to  the  piecewise  discrete  harmonic  functions  given  by 

(2.11)  a(uk,VH)  =  0,         VhGH^{n,)nV\  i  =  l,...,N 
or  alternatively  by 

(2.12)  A'l^<''  +  A-12x^'=0 

It  is  convenient  to  rewrite  (2.10)  in  variational  form.  Let  5,(u,  v)  axid  s(u,  v)  denote  the 
forms  defined  by  (2.9),  i.e. 

(2.13)  ^,(u,iO  =  a:£^^''\?      and      s(u,v)  ^  xlSys 
Equation  (2.10)  can  then  be  rewritten  as 

(2.14)  s{u,v)  =  (g,v)L2^r)  ,      veV^iT) 

Here  u  is  the  discrete  harmonic  part  of  the  solution  and  V'^{T)  C  H^^^iT)  the  restriction 

of  ^''(fi)  to  r. 

Problem  (2.14)  can  be  solved  by  different  iterative  methods  of  additive  Schwarz  type. 

3.  An  abstract  additive  Schwarz  method.  We  now  describe  and  analyze  the 
convergence  of  abstract  additive  Schwarz  methods.  We  note  that  the  theory  presented 
here  is  a  modified  version  of  the  theory  developed  in  our  previous  papers,  see  [5],  [6],  [7]. 
The  difference  is  that  we  now  include  the  effects  of  inexact  solvers  from  the  very  beginning. 

Let  V  be  a  finite  dimensional  space  with  the  scalar  product  a{u,v).  We  consider  the 
abstract  problem 

(3.1)  a{u,v)  =  f(v),  veV 


Let 

V  =  Vo  +  Vx  +  ■  •  •  +  T'V      ■ 

and  let  the  bi{u,v)  be  symmetric,  positive  definite  bilinear  forms  on  V,  x  V,.  We  introduce 
projections  T,:  V  — +  Vi  by 

(3.2)  kiTiU,v)  =  a{u,v)  .  v  E  Vi 

and  put 


We  replace  (3.1)  by 


T  =  To  +  Ti  +  •  ■  ■  +  T/v 


:=0 

Theorem  3.1.  Let 

(i)  there  exist  a  constant  Co  such  that  for  all  u  G  V  there  exists  a  decomposition 
u  =  ^i=oU,,  u,  G  K,  such  that 

N 
2_J  &i(iii,  Ui)  <  Coa(u,u) 

1=0 

(ii)  there  exist  a  constant  uj  such  that  for  i  =  0, .  . . ,  N, 

a(u,ii)  <  LL>bi{u,u)  ,  u  E  Vi 

(iii)  there  exist  constants  e,j,  for  i,j  =  1,.  . .  ,N,  such  that 

a{ui,Uj)  <  eija{u,,u,)^^^ a{Uj,Uj)^^^  ,      u^  6  Vj  ,      Uj  G  1^ 

Then 

CQ^a{u,u)  <  a{Tu,u)  <  (p{e)  +  l)uj  a(u,u)  ,  u  G  V 

where  p{e)  is  the  spectral  radius  of  the  matrix  e  =  {^o}i\=i- 
A  proof  of  this  theorem  can  be  found  in  [8]. 

4.  A  wirebasket  based  method.  In  this  section,  we  describe  and  analyze  an  iterative 
substructuring  method,  recently  developed  by  B.  Smith,  see  [9],  [10].  For  the  description 
and  analysis  we  use  the  general  framework  of  Section  3.  Here  V  =  V'^(T)  and 

(4.1)  a{u,  v)  =  s{u,v) 

4.1.  A  method  without  a  vertex  space.  To  describe  the  method,  we  need  some 
notations.  Let  T,j,  be  the  open  faces  and  let  W,  be  the  Avirebasket  of  the  substructure  f2,. 
jF,j  is  the  closure  of  the  face  common  to  the  substructures  f2,  and  Q,j.  The  wirebasket  is 
the  union  of  the  closures  of  the  edges  of  f],.  Let  Wth  and  d!Fijh  denote  the  set  of  nodal 
points  belonging  to  W,  and  dTtj,  the  boundary  of  J^ij,  respectively.  Let 

(4.2)  ui  =  —    Yl    ^(^)  '  "'J  =  —     H     "^^) 


where  n,  and  n,j  are  the  number  of  nodal  points  on  W,  and  5^,j,  respectively.  We  can 
now  introduce  the  representation 

(4.3)  V'{T)^Vo{T)  +  J2V,,{r) 

Here  VijiT)  is  the  space  of  functions  v  G  F''(r)  which  vanish  at  all  nodal  points  not  on 
J^ij.  The  space  Vq  is  defined  in  a  special  way.  It  is  a  space  of  continuous,  piecewise  linear 
functions  defined,  on  5Q,,  by  its  values  given  on  the  wirebasket  W,  and  by  constant  values 
at  the  nodal  points  of  each  face  ^,j.  The  constant  associated  with  ^,j  is  the  average  of  the 
nodal  values  on  dJ^ij .  The  basis  functions  of  Vq  are  of  the  form 

(4.4)  (^|''(x)  =  <^fc(.r)  +  — e„(.r)      on      J^„ 

Here  4>k{^)  is  the  standard  nodal  basis  function  for  Xk  G  W,  and  ^,j(x)  G  V^  is  given  by 

„  J  1      at  the  nodal  points  x  G  T,j 

'■'  I  0      at  the  nodal  points  x  G  dJ^,j 

It  is  also  possible  to  consider  the  coarse  space  Vq  as  the  range  of  an  interpolation  operator 
defined  on  Tij  by 

(4.5)  IhU=      ^      u(xk)(pk{x)  +  Ui.ie,j{x) 

We  now  define  the  quadratic  forms  corresponding  to  the  different  subspaces.  Let  bo{u,  v):  Vq 
X  Vo  ^  i?  be  of  the  form 

(4.6a)  bo{u,v)  =  {l  +  \og-)J2    E    h{u(x)-u,){v{x)-v,) 


•     X6W.,, 


and 


(4.66)  bo{u,v)^(l+\og^)Y,P,    E    hiuix)-u,)iv{x)-v,) 

for  model  problems  I  and  II,  respectively.    To  define  6,j(u,f):  Vij  x  Vij  — ►  R,  let  Q,j  = 
Q,,  U  Qj  U  J-'ij  and  introduce  the  discrete  harmonic  extension  operator  ?f,j  from  J^,j  to  J),j 

by 

.4  ^.  au.An^jU.v)  =  0 ,       ve  F''(a) u  V'in,) 

TiijU    =  u      on   J^,j  ,  HijU  =  0  on  dCl,j 
Here  V"'(Q.)  =  ^''(fi)  n  H^i^,)  and  ?^,jU  G  V'''(f^,j).  We  define  the  bihnear  form  as 

(4.8)  b,j{u.v)^  a{H,jU,n,jV) 


We  have  now  defined  all  our  subspaces  and  the  associated  bilinear  forms.  As  in  the  abstract 
theory,  we  define  the  projections  by 

(4.9a)  boiToU,!')  =  s{u,v)  ,  v  G  Fo 

and 

(4.96)  bij{T,jU,v)  =  s(u,v)  ,  v  e  V,^ 

and  let 

T  =  To  +  ^  T, , 

'J 

Theorem  4.1.  For  all  u  G  V'''(r),  with  To  and  T,j  defined  by  (4.9),  we  have 

jj 

(4.10)  7o(l  +  \og—)-K^{u,u)  <  s{Tu,u)  <  j,  s{u,u) 

Here  70  and  71  are  constants  independent  of  H  and  h  and  the  jumps  of  the  coefficients. 

To  prove  this  theorem,  we  use  Theorem  3.1  of  Section  3;  cf.  Dryja,  Smith  and  Widlund 
[4].  For  an  alternative  proof,  see  Smith  [9],  [10]. 

We  now  briefly  discuss  how  the  method  can  be  implemented.  Problem  (2.14)  has  been 
replaced  by 

(4.11)  Tu  =  g 

where  g  =  go  +  ES'ij.  9o  =  TqW/i,  and  gij  =  T./u/i- 

To  solve  (4.11)  we  can  use  the  conjugate  gradient  method  since  T  is  symmetric,  positive 
definite  and  well  conditioned.  For  simplicity  we  only  consider  the  first  Richardson  method. 
Thus  ^ 

^,"+1   =  i,n  _  ^rp^^n  _  y)   ^  ^  ^  2/(70(1  +  log  -r)-'  +  7l) 

n 
Let  r"  =  r(u"  -  tz)  =  rS  +  Eij  rf^,  where  r^  =  ro(u"  -  u)  and  rf^  =  r.j(u"  -  u).  To  find 
Tq  and  r^j,  we  solve 

(4.12a)  boir^,v)  =  ^(u", u)  -{g,v)  =  F{v)  ,  v  e  Vo 

and 

(4.126)  b,j{r^^,v)  =  Fiv)  ,  vGV„ 

To  compute  s{u"-,v)  we  solve  the  Dirichlet  problems 

(4.13)  aiHkU^v)  =  0  ,  ve  H^iflk)  n  V''{Q) 

with  li"  given  on  5fifc.  Htu"  is  the  discrete  harmonic  extension  of  u"  from  ^^a:  to  Qk-  The 
problem  (4.12a)  reduces  to  a  system  with  a  sparse  matrix  and  block  Gaussian  ehmination 


is  used  to  find  r^ix)  at  the  nodal  points  of  x  6  W,;  cf.  Dryja,  Smith  and  Widlund  [4].  To 
find  r,",  we  solve  a  problem  similar  to  (4.7).  after  replacing  Qk  with  fi.j. 

4.2.  The  method  with  vertex  spaces.  In  this  subsection,  we  discuss  a  variant  of  a 
method  discussed  in  Chapter  4  of  Smith  [9].  The  estimate  in  Theorem  4.1  contains  two  log 
factors.  We  can  remove  one  of  them  for  model  problem  I  by  adding  vertex  spaces  to  the 
representation  (4.3)  of  V^{T);  we  call  this  a  wirebasket  based  method  with  vertex  spaces. 
For  each  substructure  vertex  Xk,  we  define  a  space  VkiT)  as  follows.  Let  j^^  be  the  union 
of  parts  of  the  faces  which  have  the  vertex  Xk  in  common.  We  assume  that  dist(xfc,  dTk)  is 
on  the  order  of  H.  The  space  VkiT)  is  the  subspace  of  functions  belonging  of  V"''(r)  which 
vanish  at  the  nodal  points  not  in  J-k-  Clearly 

Let  TikU  be  the  extension  of  u  from  ^k  to  Clk  here  Qk  is  a  neighborhood  of  Xk  which  which 
contains  J^k-  Let 

(4.14)  bkiu.v)  =  a{HkU,Hkv)  ,  u,y  e  ^(r) 
where  a{u,v)  is  given  by  (2.2).  Let  Tk  be  given  by 

(4.15)  bk{TkU,v)  =  s{u,v)  ,         veVk{T) 
and  let 

(4.16)  r  =  To  +  ^r,  +  x:T., 

k  ij 

Theorem  4.2.  For  all  u  G  V'iT),  with  To,  T.j  and  Tk  defined  by  (4.9)  and  (4.15)  and 
for  model  problem  I, 

jj 

7o(l  +  log  -T-)~^  s{u,  u)  <  ${Tu,  u)  <  7i5(u,  u) 
h 

where  70  and  ji  are  constants  independent  of  H,  h. 

A  proof  of  this  theorem  is  based  on  Theorem  3.1  of  Section  3;  cf.  Dryja,  Smith  and 
Widlund  [4]. 

5.  The  Neumann-Neumann  method.  In  this  section,  we  give  a  description  of  a 
method,  introduced  by  R.  Glowinski,  P.  Le  Tallec,  Y.-H.  De  Roeck  et  al.,  see  [1],  [3],  as  an 
additive  Schwarz  method.  We  then  extend  the  method  to  the  case  of  a  large  number  of 
substructures  introducing  a  coarse  space  Vq,  which  is  similar  but  not  identical  to  that  of 
Section  4.  We  also  modify  this  method  to  handle  problems  with  discontinuous  coefficients. 

5.1.  The  method  writhout  a  coarse  space.  We  first  consider  model  problem  I, 
i.e.,  when  s{u,v)  corresponds  to  the  bilinear  form  a{u,v)  given  by  (2.2).  Let  V^(r)  be  the 


subspace  of  functions  v  6  V^iT)  which  vanish  at  all  nodal  points  on  F  \  50,.  It  is  easy  to 
verify  that 

(5.1)  y''(r)  =  i'i(r)  +  ---  +  vv(r) 

Let  Si{u,v)  be  the  bihnear  form  on  V,  x  Vi  defined  in  Section  2, 

(5.2)  .,(u,r)  =  -r^>^S<Si' 

For  the  method  considered,  using  the  notations  from  Section  3,  V  =  V'^(T),  Vi  =  K(r), 
i  =  I, . . .  ,N,  and  a{u,i->)  =  s(u,v).  There  is  no  coarse  space  Vq.  To  define  6,(u,i'):  K(r)  x 
ViiT)  — >  R,  we  introduce  a  counting  function  a,  G  V'^{T),  associated  with  9f2,. 

[2        .T  G  J-,,;, 

Q,(X)   =    <     Tltix)      X  e  W'i.ft 

[       0         all  other  nodal  points  on  F  . 

Here  n,(x)  is  the  number  of  substructures  which  have  the  nodal  point  x  6  W,  in  common. 
Let 

(5.3)  bi{u,v)  =  s,ilh(a,u)  ,  h{a,v) )  ,      u,v  e  V'iT) 

Here  //,  is  the  linear  interpolation  operator  on  the  fine  triangulation.  We  note  that  bi{u,v) 
is  symmetric  and  positive  definite  on  Vj  x  Vi. 
We  introduce  a  projection  T,:  V^  —*  Vi  by 

(5.4)  bi{T,u,v)  =  s{u,v)  ,         veVi 

and  put 

T  =  Ti  +  •  • .  +  Tat  . 

Theorem  5.1.  For  all  u  €  V^{T),  with  T,  defined  by  (5.4)  and  for  model  problem  I 


(5.5)  losiu.u)  <  s{Tu,u)  <  71  1^2  ^("^") 


(l  +  lo£if//,)2 
i/2 


where  70  ariti  71  are  constants  independent  of  H  and  h. 

This  result  is  proved  in  [3]  using  other  tools.  We  believe  that  our  proof,  which  is  based 
on  Theorem  3.1  from  Section  3,  is  simpler;  cf.  [8].  We  note  that  an  estimate  such  as  (5.5) 
but  with  four  log  factors  was  given  already  in  [4].  For  the  case  when  there  is  a  red-black 
ordering  of  the  substructures,  we  also  showed  that  three  log  factors  are  enough  and  we 
also  derived  similar  estimates  for  algorithms  with  coarse  spaces. 

5.2.  A  method  w^ith  a  coarse  space.  We  now  describe  the  Neumann-Neumann 
method  with  a  coarse  space  Vq.  The  factor  1/H^  can  now  be  removed  from  the  estimate 
(5.5).  A  function  u  G  Vo  is  a  continuous,  piecewise  linear  function  defined  by  the  nodal 
values  on  the  wirebasket  W,  and  it  is  constant  on  the  faces  !Fi-j  with  different  constants  on 
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different  faces.  This  space  is  different  from  the  coarse  space  used  in  Section  4.  The  basis 
functions  of  Vq  are  of  the  following  form.  For  a  nodal  point  Xf;  G  W,,/j,  we  use  a  standard 
nodal  basis  function  <t>k{3:)  and  for  a  face  jr,j  the  function  9,j(x)  defined  in  Section  4.  Vq 
can  be  regarded  as  the  range  of  an  interpolation  operator  defined  on  dQ,t  by 

(5.6)  hu=     Yl     u{^k)<Pk(x)+     Yl     u{i-,j)0,:(x) 

on  dCti.  Here  Xij  is  a  fixed  nodal  point  on  ^,j,  and  the  summation  in  the  second  term  is 
taken  over  the  faces  J^,j.  The  spaces  VJCF),  7  =  1 ,N,  axe  defined  by 

(5.7)  V,  =  {ve  y''(r):  v{x)  =  0 ,  X  eT\dn, ,  vf  =  0} 

Here  uf  is  a  weighted  discrete  average  of  u  on  Wj/j,  i.e. 

«r  =  —    Y    a,(x)ti(x) 
"'  xew,h 

It  is  easy  to  verify  that 

(5.8)  V''{T)  =  Vo  +  Vi  +  ---  +  Viv 

We  now  introduce  bilineax  forms  b,(u,v):  V,  x  Vi  — >  R.  For  i  =  1, . . .  ,N ,  they  axe  defined 
as  in  (5.4).  For  i  =  0,  we  use 

/ggs  bo{u,v)     =  (1  +  logf  )-iE.{/^Ex6W,h(^^(^) -"«)(^(^) -'^t) 

+  E;r,^can.  H{u{xij)  -  Ui){v{x,j)  -  u.)} 

Here  u,  is  the  discrete  average  of  u  on  Wj/i  introduced  in  (4.2).  Let 

r  =  To  +  Ti  +  •  •  •  +  T;v 

where 

(5.10)  hi{TiU,v)  =  s{u,v)  ,         u  e  V;  ,  i  =  0,...,iV 

Theorem  5.2.  For  all  u  G  F''(F),  with  Ti  defined  by  (5.10)  and  for  model  problem  I 

(5.11)  los{u,u)  <  siTu,u)  <  7](1  +\og—Ys{u,u) 

h 

where  70  and  71  are  constants  independent  of  H  and  h. 
A  proof  of  this  theorem  is  given  in  [S]. 

5.3.  A  method  with  a  coarse  space  for  model  problem  II.  In  this  subsection, 
we  describe  a  Neumann-Neumann  method  for  model  problem  II.  The  bilinear  form  is  now 
given  by 

(5.12)  a{u,v)^   I  p[x)Vu-Vvdx 

Ju 
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where  p(x)  =  p,  =  constant  >  0  on  substructure  Q,.  Let  Vi,  i  =  1,. . .  ,N,  associated  with 
dClt,  be  a  space  of  functions  i'  £  F''(r)  which  vanish  on  the  wirebasket  W,  of  fi,  and  at  all 
nodal  points  on  F  which  are  not  on  dO,,.  The  space  Vq  is  the  same  as  in  Section  4  i.e.  it  is 
the  range  of  the  interpolation  operator  Ih]  ci.  (4.5).  We  have 

(5.13)  v''ir)^Vo  +  \\  + ■■■  +  ¥!, 

To  define  b,{u,v):  V,  x  V,  —*■  R,  i  =  1 N,  we  introduce  a  piecewise  constant  function 

Pi,  defined  on  dQi,  which  is  equal  to  {p,  +  Pj)^^^  on  the  face  Tij  common  to  f2,-  and  Clj.  Let 

(5.14)  k{u,v)  =  s,(h{p,u)Jh{piu) )  ,      u.v  E  V, 

where  Si{u,v)  corresponds  to  the  form  a{u,v)  —  ( Vu,  Vu)z,2(n);  cf.  (2.2)  and  (2.13).  For 
?  =  0 

(5.15)  bo{u,v)  =  Y^p,    J2    h(u(-r)-Ui)(v(x)  -V,) 

where  u,  is  defined  in  (4.2).  Note  that  for  u  G  Vj,  Ik{p,u)  is  defined  uniquely  on  the 
wirebasket  W,  since  u  =  0  there.  Let  T,:  V^{T}  — >  Vi  be  a  projection  defined  by 

(5.16)  b,{TiU,v)^s(u,v)  ,         V  6  K  ,      i  =  0,l,...,iV 

and  let 

T  =  To  +  Ti  +  •  •  •  +  Tn 

Theorem  5.3.  For  all  u  G  V'''(F),  with  the  Ti  defined  6j/ (5.16),  and  for  model  problem 
II 

H  H 

(5.17)  -foil +\og— )-^s{u,u)  <s{Tu,u)  < -^{1  +  log— fs{u,u)  . 

n  h 

Here  70  and  71  are  constants  independent  of  H  and  h  and  the  jumps  of  p{x). 
A  proof  of  this  theorem  is  given  in  [S]. 

6.  A  multilevel  method.  In  this  section  we  discuss  a  modification  of  the  multilevel 
method  suggested  in  J.  Bramble  et  al.  [2]  for  the  problem  (2.14)  in  terms  of  the  additive 
Schwarz  method  described  in  Section  3.  We  first  consider  model  problem  I  and  then 
model  problem  IL  We  note  that  we  have  previously  described  this  multilevel  method  as 
aji  additive  Schwarz  method  for  the  original  problem  (2.3);  cf.  Dryja  and  Widlund  [7]. 

6.1.  The  method  for  model  problem  I.  We  discuss  problem  (2.3)  with  the  bihnear 
form  given  by  (2.2)  defined  on  a  triangulation  of  f2  obtained  by  successive  refinements.  We 
consider  (£  +  1)  levels  of  triangulations  of  Q.,  associated  with  parameters  hk,  which  satisfy 
y''*-'(f))  C  V^^{9,).  Here  the  V^>'{?L)  are  standard  finite  element  spaces  with  ¥^"{^1)  and 
V^'[Q,)  the  coarsest  and  finest  spaces,  respectively.  V^^'iT)  is  the  restriction  of  V^''{Q.)  to 
F.  The  problem  considered,  (2.3),  is  defined  in  F'"(F).  As  a  representation  of  V''  =  V''', 
we  take 

(6.1)  V''  =  Vo  +  V,  + ■■■  +  ¥,■  ,         Vu  =  V^' 

12 


We  introduce,  for  k  >  1,  bkiu.v):  Vk  x  T4  — +  R  by 

(6.2)  bkiu,v)=   Yl  hku(x)v(x) 

xeNk 

Here  A'^^  is  the  set  of  nodal  points  of  the  level  A-  triangulation  on  T.  We  introduce  T^:  V'^  —y 
nby 

(6.3a)  s{Tou,v)  =  siu,v)         v  e  Vo  , 

(6.36)  bkiTkU,v)  =  s{u,v)  ,  v  e  Vk  ,  k  =  l,...,N 

and  put 

r  =  To  +  Ti  +  •  •  •  +  Tf 

Theorem  6.1.    Let  hk/hk+i  be  uniformly  bounded.     Then  for  all  u  6  ^''(r),  with  Ti 
defined  by  (6.3),  and  for  model  -problem  I 

(6.4)  7o(^+l)"^s(u,u)  <  s[Tu.u)  <  fi(e  +  l)s(u,u) 
If  Cl  is  convex  then 

(6.5)  725(u,u)  <  s{Tu,u)  <73(^+l)s(u,u) 

Here  7,,  i  =  0, .  . . ,  3,  are  constants  independent  of  Hq,  . . .  ,h(. 

Proof:  The  proof  uses  Theorem  3.1  of  Section  3. 

Assumption  (i):  Let  Hu  be  the  discrete  harmonic  extension  in  the  sense  of  a(u,  v)  of 
u  from  dQi  to  n,,  for  all  i.  We  use  the  partitioning 

e 
u  =  J2^'=^      uo\r  =  Qonulr  e  Vo(r)  ,      Uk\r  =  (Qk  -  Qk-i)nu\^  e  Vk{T)  ,  k  =  !,...,£ 

k=0 

Here  Qk  is  the  L2-projection  from  V'^(Q.)  — >  Vk{Q),  i.e., 

iQkU,v)L2(Q)  =  {u,v)mu)  ,     veVk,  k  =  o,...,e 
It  is  well  known  that 

(6.6)  |Q.w|//i(n)  <  C\u\hi(u) 
and 

(6.7)  ||u  -  Q,u\\L^a)  <  hi\u\H^n) 
Using  (6.6),  we  have 
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For  k  =  1 £,  we  obtain,  by  using  an  elementary  inequality  and  (6.7), 

<  Ch^HWHu  -  Q^:■Hu\\l,^^^  +  \\Hu  -  Q,_i?^u||^,,„,} 

<  C\nu\jji^^)  <  Cs{u,u) 

Thus 

N 
5(Uo,Uo)  +  I^i.(tinli.)   <   C{e+l)s(u,u) 
1=1 

which  proves  Assumption  (i). 

For  a  convex  region  fl,  we  can  instead  use  uq  =  PqHu  and  Uk  =  {Pk  —  Pk-i)'Hu,  where 
Pk  is  the  iJ^ -projection  and  use  the  fact  that  the  resulting  Uk  are  /f^ -orthogonal;  cf.  Dryja 
and  Widlund  [5].  For  more  details  see  below. 

Assumption  (ii)  follows  from  the  inverse  inequality. 

Assumption  (iii)  is  proved  by  using  standard  arguments. 

6.2.  A  new  method  for  model  problem  II.  We  now  discuss  a  new  method  for 
model  problem  II,  i.e.  when  the  bilinear  fox'm  is  given  by 

a{u,v)  =   /  p{x)Wu  ■ '^v dx  , 

where  p[x)  =  pi  —  constant  >  0  on  Q,,.  The  decomposition  of  V'''(r)  is  now  of  the  form 
(6.8)  y^F)  =  Vo(r)  +  V,{T)  +  •  ■  •  +  Vt{T) 

Here  Vq  is  the  space  from  Section  4,  i.e.  the  range  of  the  interpolation  operator  Ih\  cf.  (4.5). 

We  let  Vk,  /:  =  1,  ...,£,  be  the  space  of  functions  in  V^''{T)  which  vanish  at  the  nodal 
points  of  the  wirebasket  of  5f2,.    It  is  easy  to  verify  that  the  representation  (6.8)  of  V 
holds. 

We  introduce  bk{u,v):  Vk{T)  x  ^^.(F)  ->  i?,  by 

jj 

(6.10)  bo{u,v)  =  (l  +  \og—)Y,     ^    /ip,(u(a:)-u.)(u(ar)-v.)  , 

cf.  (4.6),  and  for  k  =  !,...,£, 

(6.11)  6/,(u,i;)=    X^  /ifc/5(x)K(.T)j;(a,-)  ,      u,v  eVk{T) 

x^Nk 

Here  p{x)  is  piecewise  constant  and  equal  to  (p,  -|-  /9j)/2  on  the  face  J^i^  which  is  common 
to  the  substructures  Q,  and  9.y  Note  that  for  u  6  V^,  p{x)u[x)  =  0  on  W,  since  u{x)  =  0 
there. 

We  now  introduce  T/,:  ^''(F)  ->  Ft,  by 

(6.12)  hk{TkU,v)  =  s{u.v)  ,         v^Vk  ,  k  =  Q,...,(. 
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and  put 

(6.13)  r  =  ro  +  Ti  +  ---  +  T^ 

Theorem  6.2.    Let  hk/hk+\   be  uniformly  bounded.     Then  for  all  u  G  F''(r),  with  Ti 
defined  by  (6.12),  and  for  model  ■problem  II 

(6.14)  7o(C+l)"\s(if,u)  <  s{Tu,u)  <-ti{i  +  l)s{u,u) 

where  70  and  71  are  constants  independent  of  ho,  hi, . . .  ,h(  and  the  jumps  of  p{x). 

Proof:  We  again  use  Theorem  3.1. 

Assumption  (i):    Let  uq  =  IhU  where  I^u  is  defined  in  (4.5).    It  can  be  estabhshed 
that 

(6.15)  bo(uo,Uo)  <  C{C  +  lfs{u,u) 

see  Dryja,  Smith  and  Widlund  [4]  and  Smith  [9].  We  now  define  u^.  Let  w  =  u  —  uq  and 
note  that  it  vanishes  on  the  wirebasket  W,.  The  function  w,  given  on  J-ij,  the  face  common 
to  n,  and  rij,  is  extended  to  f2,  as  a  discrete  harmonic  function  Uij,  i.e., 

(6.16)  (Vii;,,,  Vi;)z.2(o.)  =  0  ,  ve  V^{n,) 

Wij  =  w   on   Tij    and   uuj  =  0   on   dO.,  \  J-ij 

Let  V'"'{Q.,),k  =  0, ■■-,£,  be  the  subspace  of  V^^M^i),  the  restriction  of  V^"  to  fl,,  of 
functions  which  vanish  on  dQi\Tij.  We  note  that  il\j  6  V'^'{Qi).  This  function  is  represented 
as 

(6-17)  Wij  =  Pij,iw,j  +  (P,j,2  -  P,j,i  )w^J  +  •  •  •  +  {Pij^e  -  P,j,e-i  )wtj 

where  P.j,^  is  the  if  ^-projection:  V'^i^,)  -^  V'"'{Q,i),  i.e. 

(VP,,,,v,V<|>)L.^a.)  =  (Vi;,  V<^)i2(f,.)  ,      <(>  e  V'H^,) 

It  is  known  that,  since  fi,  is  convex, 

(6.18)  \\u  -  P,_,,fct<||£,2,n.)  <  Mhk\u\fmQ,) 

On  J^tj  we  define  u^  as  the  traces  of  the  functions  defined  in  (6.17) 

ui  -  P,j,\Wij  ,      Uk  =  (P,j,k  -  Pij.k-i  )uuj  ,  k  =  2 , £ 

It  is  easy  to  see  that  P,j,yt  —  Ptj,k-\  is  an  orthogonal  projection  from  V'^{Q,i)  into  l''''*(fi,)  = 
Ra.nge{Pij,k  -  Pij,k-i)  C  V'"'{fl,).  We  now  show  that 

e 

(6.19)  E^^-("'t,iu)<C'(6'  +  l)'5(u,u) 
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for  the  Uk  defined  above.  We  obtain,  by  using  the  definition  of  p{x), 

(6.20)  b,{u,,u,)=   Y,  h>.p{x)icl{.v)^lY.hkJ2     E    P'^li^) 

where  Nij^k  is  the  set  of  nodal  points  of  the  A:-level  triangulation  of  JF^  and  the  summation 
with  respect  to  j  is  over  the  different  faces  of  dQi.  For  one  face,  we  obtain 

(6.21)  <  Chfpi{\\P,j,,u,  -  u,||^,(„,)  +  \\P,j,k-iUk  -  tU-||L(n.)} 

<  Cp.luklHuu.) 

Adding  these  inequahties  with  respect  to  k,  we  obtain 

PiEi=i  hkY,TeN„,kUl{x)  <  C p^Yli=i{^Uk,'^Uk)mQ,) 

=  C'p.{(VP,j,iu;,j,  Vu;,_,)i,2(n,)  +  T,i=2i^iP'j,k  -  P,j,A:-i)u;.j,  Vu;,j)i,2(n,)} 
=  Cpi{\'Pij,eUij,Vu!ij)mu,)  =  Cp,\\Vu!ij\\l^^^.^ 
<Cp.iu.'|2   /  <C(i+l)'s,{u,u) 

For  the  last  inequality,  see  for  example  Dryja  and  Widlund  [6].  Combining  this  inequality 

with  (6.20),  we  obtain 

(. 

k=\ 

Together  with  (6.15),  this  provides  the  bound  for  Assumption  (i). 
Assumption  (ii):  It  is  known  that  for  u  G  Vo 

a{u,u)  <  Cbo{u,u) 

cf.  [4].  For  k  =  l,...J, 

a{u,u)  <  Cbi;{u,u)  ,  u  £  Vk 

follows  from  an  inverse  inequality. 

Assumption  (iii)  is  proved  by  using  standard  arguments. 
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